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EQUIVALENCE OF WEAK GALERKIN METHODS AND VIRTUAL 
ELEMENT METHODS FOR ELLIPTIC EQUATIONS 


LONG CHEN 


Abstract. We propose a modification of the weak Galerkin methods and show its equiv¬ 
alence to a new version of virtual element methods. We also show the original weak 
Galerkin method is equivalent to the non-conforming virtual element method. As a conse¬ 
quence, ideas and techniques used for one method can be transferred to another. The key 
of the connection is the degree of freedoms. 


1. Introduction 

Recently several numerical methods |E 13 0 [2| have been developed for polygon and 
polyhedral meshes. In this paper, we shall discuss the connection between weak Galerkin 
(WG) methods and virtual element methods (VEM). We show that, for diffusion equations, 
a modified WG can be derived from a new version of VEM or equivalently a new version 
of VEM can be derived from WG. In the same sprit, we will show the original version of 
WG J7) is equivalent to non-conforming VEM j4). The equivalence enables us to apply the 
convergence theory, as well as computer implementation, of VEM to WG and vice verse. 
It should also help in giving new insights for each method when applied to other equations. 

The key of the connection is the degree of freedoms (d.o.f.). VEM space can be em¬ 
bedded into WG space through the degree of freedoms. Actually WG space contains more 
d.o.f. but some can be eliminated locally and contribute to the stabilization term only. Es¬ 
sentially WG and VEM share the same degree of freedom but the associated functions are 
different. In VEM, the shape function is determined by a suitable PDE inside each element 
which is in general non-polynomial. The point-wise information of the shape function, 
however, is not needed in the computation which leads to the name ‘virtual element’ . The 
resulting space is conforming. In WG, inside each element, the shape function is simply 
polynomial. These functions are totally discontinuous across elements. The continuity is 
imposed by the stabilization through a suitable boundary integral defined on the boundary 
of elements. Namely, in WG, we know more on the shape function by sacrificing the con¬ 
tinuity. The violation of the conformity is under control as we will show WG will also pass 
the patch test and recast in terms of VEM. 

The rest of this paper is organized as follows. We first present the two methods follow¬ 
ing the presentation in the literature. We then give a modified form for each of the methods 
and then show they are identical. We will try to stick to the notation used in each method 
such that the readers who are familiar with one or another can easily follow. In particular, 
we use letter V for spaces in VEM and W for those in WG. 


Date : March 17, 2015. 

LC has been supported by NSF Grant DMS-1418934. 


1 



2 


LONG CHEN 


2. Function Spaces 

We consider a two dimensional domain 12 which is decomposed into a polygon mesh 
Th- Each element is a simple polygon and denoted by K. We use two dimensional case for 
the clear illustration and will talk about the generalization to high dimensions afterwards. 

2.1. Spaces in Virtual Element Methods. We introduce the following space on K 

(1) V k (K) := {u 6 H\K) : v\ aK G M k (dK),Av G P k - 2 (K)}, 

where IP/, ( D ) is the space of polynomials of degree < k on D and conventionally P_j ( D) = 
0, and the boundary space 

B k (dK) := {v G C°(dK) : v\ e G Pfe(e) for all e C dK}. 

The shape function is defined in ([[} but the point-wise value of a function v G V k (K) 
requires solving a PDE inside K and thus considered as implicitly defined not explicitly 
known. The novelty of VEM is that not the point-wise value but only the d.o.f. is enough 
to produce an accurate and stable numerical method. 

Consider the dual space 

(2) span {x v ,Xe~ 2 >X k K 2 }’ 

where the functional vectors 

• x v : values of v at the vertices of K ; 

• Xe~ 2 '- the moments on edges 

Xe{v) = \e\~ 1 (m,v) e Vm G M fc _ 2 (e), V edge e C dK; 

• Xk : the moments on element 

Xk(v) = \K\~ 1 (m,v) K Vto G Mfc_ 2 (iv). 

Here the scaled monomial M/( D) on a (/-dimensional domain D is defined as 

with x c the centroid of D and hjj the diameter of D or hp = |Z?|t/ 2 . 

The verification 

(3) (V k (K))' = s P an{ Xl „ X ^ 2 }, 

is called unisovlence and can be found in |2j. 

Remark 2.1. The PDE operator A used in the definition of VEM space 0} can be re¬ 
placed by other operators as long as f‘ k (K) C V k (K), which ensures the approximation 
property of V k {K), and the unisolvence holds. For example, when K is triangulated into a 
triangulation Tk , we can chose the standard k -th order Lagrange space on Tk and impose 
A hV G Pfc_ 2 (A') where A/, is the standard Galerkin discretization of A. 

We relabel the d.o.f. by a single index. Associated to each d.o.f., there exists a basis of 
V k (K) such that Xii^j) = $ij f° r *> j = 1,..., N k = dimf4(A'). Then every function 
v G 14 (K) can be expanded as 

N k 

V = ^2xi(v)(l>i 

i= 1 

and in numerical computation it can be identified to the vector v G M. Nk 

v = {Xi(v),X2(v), ■ ■ ■ ,XN k {v)). 
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The isomorphism can be denoted by 

X : V k (K) ->• R Nk , x{v) = (xi(«).X2(«), ■ ■ -,XN k (v)). 

The inverse of this isomorphism will be denoted by 

<1> : R Nk —»• Vk(K), $(v)=4>-v, 

if we treat the basis = {cf> i, fa,., (f>M k } as a vector. 

Later on, to compute the L 2 -projection of the shape function, the authors of fT] intro¬ 
duce a larger space 

(4) Vk(K) = {w G H\K) :v\ dK & M k (dK), Av G P k (K)} 
and a subspace isomporphism to V k (K) 

(5) V*(K) = {w£V k (K) : {w-I^w,q*) K = 0 Vq* G M k (K)\M k _ 2 (K)}, 

where the projection I [J 7 will be defined in the next section. 

The dual space (V k (K))' is enlarged from (' V k (K ))' by including volume moments of 
order k — 1 and k. The spaces V k (K) and (K) are different but share the same d.o.f. 
For the same vector v G K' Vfc , we can then have different functions associated to them and 
in general ^ c l>yR(ii). 

Function spaces in each element will be used to design a virtual element space on the 
whole domain. Given a polygon mesh 7/t of O and a given integer k > 1, we define 

V h = {vsH 1 (Q):v\ K €V k (K) VKcT h }, 

V h = {vsH 1 (Q):v\ K €V k (K) VK c%}, 

V? = {v G H\n) : v\ K G V^(K) MK C %}. 

The d.o.f. can be defined for the global space in the natural way. 

For pure diffusion problem, V k is enough. The introduce of and Vy, will be helpful to 
deal with low oder terms, e.g., reaction-diffusion problem, and simplify the implementation 
in three dimensions. 

2.2. Spaces in Weak Galerkin Methods. We hist introduce the weak function space 

(6) W{K) ■= {v = {u 0 ,u b } : v 0 G L 2 (K),v b G L 2 (dK)}. 
and, for k > 1, the weak Galerkin finite element space 

W k (K) := {d = {uo, v b } : vq G P k {K),v b G B k (dK)}. 

We further split the space as the summation of boundary and interior part 

W k (K) = W 0 (K) + W b (K), 

where 

Wo(K) = {ve W k (K) I V = {u 0 , 0}}, 

W b (K) = {»G W k {K) | v = {0, w 6 }}. 

For any function v = {« 0 ; v b } G W k (K), it is convenient to extend the notation of v 0 and 
v b so that, without ambiguity, vq G Wq(I\ ) and v b G W b {K). 

The space H 1 ^) can be embed into the weak space W(K) and project onto W k {K) 
by the L 2 -projection as 


Q h u := {Qou , Q b u}, 
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where Qq : L 2 (K) —s- Pfc(iT) and Qb : L 2 (dK) —> B k{dK) are L 2 -projections. The 
projection Qh will preserve P k{K) in the sense that 

QhU k = {u k , ilk \dK } for all u k G P, k(K). 

We will simply write as P k{K) C Wk(K). 

The global weak Galerkin finite element, for a given polygon mesh Th and fc > 1, is 
defined as 

(7) W h := {v G W : v = {v 0 ,v b },v\ K G W k {K) VK C %}. 

Remark 2.2. In the original WG space, the boundary part is defined as Vb\ e G Pfc(e), Ve C 
dK and can be further reduced to Vb\ e G Pfc_i(e),Ve C dK. That is vi, is also discontin¬ 
uous at vertices. In <0, Uf, is continuous on the skeleton of 74. We will come back to the 
original version and show the connection with the non-conforming VEM. 

Obviously we can identify the dual space 

{Wk(K))' = (Vk(K)Y = span {x v ,Xe~ 2 ,Xk} = span{x b ,x 0 }, 
which introduces an isomorphism 

ly-j -W := ° X '■ Vh -> Wh- 

The space Vh can be embed into Vj, by the zero extension of corresponding vectors. That is 
for atG Vh, we can identify it as a function in 14 with the same d.o.f. except the volume 
momentum of order k — 1 and k are zero. Consequently Ir, yW \v h leads to an embedding 
of I v->w '■ Vi Wh and will be abbreviated as I. Without ambiguity, the isomorphism 
I w^v also abbreviated as I. 

In the implementation, the practitioners will use the same vector and the same inner 
product in Euclidean spaces. In the continuous level, however, the same vector could link 
to different functions in different spaces and thus leads to different interpretation. 

Remark 2.3. The WG space Wh corresponds to a larger VEM space V/, . The space 14 
and 14 r will be isomorphism to a reduced space to be discussed later. 

To be consistent to the WG space, we classify the d.o.f. to boundary part \b an( f interior 
part x 0 - Note that Xb( v o) = Xb({^o, 0}) = 0 but Xb(v 0 \ d K) Y 0. 


3. Methods 

We use the simplest Poisson equation for a clear illustration. The equivalence of VEM 
and WG will hold for general diffusion equations. We will also comment on their differ¬ 
ence which lies in the treatment of the low order terms. 

Consider the Poisson equation with zero Dirichlet boundary condition: 

—Am = / in rtlan = 0. 

The weak formulation is: find u G //,| (G) such that 

a(u,v ) := (Vu, Vv) = (f,v) Vu G Hg(fl). 
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3.1. Stiffness Matrix. We can chose a conforming virtual finite element space V k := 
Vi, n //(! (12). However, we cannot compute the Galerkin projection to Vf since the tradi¬ 
tional way of computing a(uh, Vh) using numerical quadrature requires point-wise infor¬ 
mation of functions and/or their gradient inside each element. In virtual element methods, 
only d.o.f is used to assemble a stiffness matrix having certain approximation property. 

Define n^ 7 : V k (K) —> P fc (AT) as the H 1 projection: given v £ V k (K), define £ 
P k(K) and X*( n ^w) = 0 such that 

(Vn^w, Vp) K = (Vv, Vp) K , for all p £ P k (K). 

The right hand side can be written as 

(Vu, Vp) K = -{v, A p) K + (v,n- '\/p) dK 

and can be computed using d.o.f. of v since, forp £ P k (K), A p £ P k _ 2 {K) and Vp • n £ 
Pfc-i(e),e £ dK. 

The stiffness matrix of the virtual element method is defined as 

a VEM (w, v) :=a(nj u ,n^). 

Remark 3.1. As Vp ■ n £ Pr- i (e), e £ dK, if we do not enforce the continuity of v\o k, 
we could discard the d.o.f. on vertices and use d.o.f. of the edge moments up to order k — 1 
which leads to the non-conforming VEM. 

Remark 3.2. The boundary integral ( v , n ■ Vp) e can be computed more conveniently if 
we use Gauss-Lobatto points as the d.o.f. on edges. 

We now turn to the weak Galerkin finite element method. Define the weak gradient 

V£ : W k {K) -A VP k {K) as: given v = {u 0 ,M £ W k {K), define £ X7P k (K) 
such that 

(8) (V^u, Vp) K = ~(vo, Ap)iu + (v b ,n ■ Vp) dK for all p £ P fc (A'). 

And the weak gradient V^, : W k —> UKc.Tj?k{K) as (V w v)\k = V^4 for all v £ Wf,. 
The stiffness matrix of WG is defined as 

a WG (u, v) := (V w u,V w v). 

Remark 3.3. The original weak gradient introduced in (TJE] is defined as V f : W k (K) —>■ 
Pfe-iW as 

(9) (V^v,q) = ~(v 0 ,V ■ q) + (v b ,n-q) dK for all q £ P^^AT), 

which requires solving a larger local problem for k > 2. The size of equation © is 
(k + 1)(A; + 2)/2 — 1 while the size of © is k[k + 1). For higher dimensions, the saving 
is more dramatic. The current modification is motivated by VEM. 

We have the following equivalence of the stiffness matrix used in WG and VEM. 

Theorem 3.4. For u £ 14, we have 

GO) vn^u = v^(iu), 

and consequently 

(11) a VEM (it,w) = a WG (lu,Iv), for all u,v £ 14. 

Proof The verification (flOb is straightforward since u and hi share the same d.o.f. □ 
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3.2. Stabilization. The stiffness matrix alone will not lead to a stable method. The lack 
of stability can be easily seen from the VEM formulation. Since IFT( K) C Vu(K) and it 
is a strict subspace except the case I\ is a triangle, we may have a VEM (u, v) = 0 when 
v £ ker(n^)/R. Namely a VEM (-, •) alone cannot define an inner product on V®. From 
the WG point of view, the weak gradient operator V,„ may have non-trivial kernel other 
than the trivial constant kernel of the gradient operator. 

A stabilization term is added to gain the coercivity. To impose the stability while 
maintain the accuracy, the following assumptions on the element-wise stabilization term 
Sk {■ , •) are imposed in VEM 0. 

• fc-consistency: for pk £ Pk{K) 

S K (pk,v) = 0, VveV h . 

• stability: 

Sk(u,u) ~ a(u,u) for all u £ (/ — 11^)14. 

We then define 

ah EM (u,v) = a(n^u,II^) + S K {u,v). 

K&T h 

The fc-consistency will imply the Patch Test, i.e., if u £ Pfc(fl), then 
a{u, v) = a^ EM (u, v), forallu£Vfc. 

The stability will imply 

a(u,u) ~ a£ EM (u,u) for all it £ V/j. 

An abstract error estimate of VEM with stabilization satisfying fc-consistency and stability 
is given in 0. So VEM is in fact a family of schemes different in the choice of stabilization 
terms. 

In the continuous level, one such choice is the scaled L 2 -inner product 

h K 2 ( u - n fcV v - v) K 

The fc-consistency is obvious as 11^ preserve fc-th polynomial. The stability is from the 
first order error estimate of u — IlJ u: 

h^\\(I -I^)u\\ K <\\V(u-U^u)\\ K . 

In the implementation, I [j 7 is realized as a matrix and the stabilization can be chosen as 

( 12 ) S x (u,v) :=x((I ■ %((/- nj». 

That is we use the l 2 product of the d.o.f. vector to approximate the L 2 -inner product. 
The scaling is build into the definition of d.o.f. through the scaling of the monomials. The 
norm equivalence of l 2 and L 2 is due to the well known fact: the mass matrix is spectral 
equivalent to its diagonal. So far (fl2l) is the stabilization used in the current VEM. 

We now discuss the stabilization of WG. For u, v £ Wk(K), we can use 

h~K ( u b ~ u o, Vb — Vo)dK- 

One particular choice is using the d.o.f associated to boundary of elements only 

(13) S dK (u,v) := Xb( u b ~ u 0 ) ■ Xbi. v b - v 0 ). 

We then define 

a^ G (u, v ) = (W w u,V w v) + Y Sok{u,v). 

KcT h 
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When u £ Pfc(iv), uo\dK = Ub . We thus also have the patch test, i.e., if u £ Pfc(H), 
then 

a(u,v) = a^ G (u,v) for all v £ i/ 1 (iv). 

The stability of WG is, however, ensured by a different mechanism: the continuity 
across the boundary of elements. 

Lemma 3.5. Ifa^ G (u, u ) = 0, then both uq and u b are globally constant and Uo\g h = u b . 

Proof aff G (u, u) = 0 implies Sok(u, u) = 0 that is uo\ok = Ub. We now prove Vtio = 
0 from the fact V w u = 0 as follows 

(V«o, Vw 0 )t = -(«o, V • Vuo)t + (uo,n ■ Vu 0 )aT 
= (m 0 , V • Vu 0 )t + (iib, n ■ Vu 0 ) dT 

= (V w u, Vu 0 ) 

= 0. 

Therefore uq is constant in I\. So is Ub. They are globally constant since uq\qk = Ub■ □ 

The constant kernel can be further eliminated by imposing the boundary condition into 
the space by defining: 

= {v £ W h : v b \dtt = 0}. 

3.3. Right hand side. In WG, it is straightforward as vq represents the function inside 
element and thus only vo contributes to the L 2 inner product. So we simply compute 

(f, v o)- 

In VEM, we have to resort to the d.o.f. For k > 2, define Qfc_ 2 as the L 2 (K )-projection 
to P k - 2 (K). Then the exact inner product (/, v) is approximated by 

(fh,v) = ( Qk- 2 f,v) K , 

KcT h 

which can be computed using d.o.f. of v £ V k (K). For k = 1, we use 

(14) (Qof,v) K , 

where v = meaiif^;,,(c)). To achieve optimal order error estimate in L 2 - norm, more 
accurate approximation of the right-hand side (RHS) is needed for k = 1,2. 

3.4. Equivalence. The VEM method is: find Uh £ V® such that 

(vnj^, vn lv h ) + Y S ^(uh, V h ) = (f h , v h ) for all v h £ V°. 

K&T h 

The modified WG method is: find Uh £ W® such that 

(15) (V w Uh, V w Vh) + Y S 9 k(u 0 - u b ,v 0 - v b ) — (f,v 0 ), for all v b £ W®. 

KcT h 

We shall propose a new version of VEM on the larger space Vh and then show that it is 
exactly WG. The new version of VEM is obtained by using a new projector II^ : V k (K) —>■ 
P k(K) in the stabilization term. 

(16) n fc : V k (K) -»• P k (K), Xk(v) = Xk^v). 

Note that in V k (K) only the momentum of oder A; — 2 is defined. To define a function in 
Pfc(if), we need to work in the larger space V k (K). Obviously Yl k p k = Pk if p k £ P k (K). 
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We propose the following stabilization 

(17) S(u, v ) := x(u - n fcw) ’ X(v - Hkv). 

Obviously S(-, •) is ^-consistent since II/, preserves the k -th polynomial. We prove the 
stability as follows 

s(n,u) ~ h~ K 2 \\{i - n fe )u || 2 = hj 2 \\(i - n k )(i - n^)u\\ 2 < ||v(/ - i$) u \\ 2 . 

The corresponding VEM method is: find u k G V® such that 

(18) (vn^fc, Vn^Vfc) + S(u h ,v h ) = (/, n k v h ) forall v h G v£. 

We now show its equivalence with WG. First of all, by Theorem l3.4l they share the same 
stiffness matrix. We then look at the stabilization term. By the definition of the projector 
life, Xk( u ~ IL-w) = 0. The stabilization (fTTb can be simplified to boundary terms only 

S{u, V ) = Xb( u - n fcu) ■ Xb( v ~ n fc«) 

Given a function u = {uo) u b} G W k (K), we use d.o.f. to identify u as a function 
lit G 14(A'). Note that Iu 0 ^ uq since Iu 0 G V k (K) is zero on dK and in general non¬ 
polynomial while uq G P k(K). What is uqI It is Ilfc(Iu). We then have the equivalence 

X b {ub - u 0 \dK ) • Xb( v b - vo\qk) = Xb( Iu - life (lit)) • Xb( Iv ~ n fc (It;)) 

Finally the right-hand side in (IT8l > is related to WG as 

(/, Ilfelu) = (f,v 0 ). 


We summarize as the following theorem. 

Theorem 3.6. The modified WG method (1151) is equivalent to the VEM method (1181 ). 


The equivalence is, however, established for Wh and ly,. We can further eliminate some 
d.o.f. in Wh to show the equivalence to ly,. To this end, for k > 2, we split the interior part 
of a weak function uo = u\ +112 such that mi G Pfc_ 2 (A') and 0 ) = Xk~ 2 { u 1 ) for 

all the volume momentum up to order k — 2. In other words m 2 G° V'j'(A), the annihilator 
of the dual space Vj'(A') and thus in the volume integral (wo, Ap) = (ui, A p) for all p G 
Pfe(Af). By the definition of the weak Gradient, we have V „,({«(,, ito}) = V w ({ub, ui}). 
Namely U 2 does not contribute to the stiffness matrix. We can then determine u-> by solving 
the equation of stabilization term only 


Xb( u b-ui\dK ~ u 2 \dK) ■ Xb( v 2 \dx) = 0 , for all u 2 G° V k {K). 


If we introduce D\ : P k -2 —> and D 2 :° Vj. '(AT) —» as the matrix repre¬ 

sentation of Xb applied to u \ and u 2 respectively, then the matrix representation of Sqk 
is 


/ / _ jDl _£> 2 \ 
l-Df DjD 1 DfD2\ 

\-D% DlDy DjD 2 J 



-d 2 ). 


After the elimination of it 2 , the stabilization matrix for the reduced space is 


{-DIP 


~PD 1 

DlPD x 



where P — I — Z? 2 (Z?JZ? 2 ) X D\ is the projection to the ker(£ff). 
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Remark 3.7. Here we eliminate u 2 to show the equivalence. In practice, we can chose 
convenient basis to expand uq (e.g. monomials) and eliminate the whole uo element-wise. 
The resulting global system will involve u b only. Of course, such reduction can be also 
applied to VEM and known as condensation in finite element methods. 

For the right hand side (/, v), for k > 2, we could write RHS of WG as (/i, f 2 ) with 
fi = (/, i = 1, 2. To achieve optimal order in // 1 type norm, we can discard f 2 and 
only keep (/i, 0). Then the reduced system will have identical RHS as VEM. 

3.5. Difference. The difference of WG and VEM lies in the treatment of low order terms. 
We illustrate it by considering the RHS. As we mentioned before, to achieve optimal order 
in H 1 norm, we only need fi and discard f 2 . To achieve optimal order in L 2 norm, we do 
need more. Let us keep both fi and /2 in WG. We then solve 

Xb(u b - mldK - u 2 \dK) -XbMdK) = / 2 , for all v 2 G° V k (K) 

to eliminate u 2 . After the elimination of u 2 , the RHS of the reduced system is modified to 
(/o, h) with f 0 = D 2 {DlD 2 )- x f 2 and f x = fi - D\D 2 (D%D 2 )~ l f 2 . In VEM, a L 2 
projection operator 1I® is introduced [Q and a perturbed RHS is computed using 1 l ( ; j which 
is different with the current formulation. 

For k = 1, the reduced system of WG is in the form 

(A + P)u = /, 

with / = l) 2 (l )2 I)>y [ f 2 . Since the optimal /A-error estimate, for k > 1, has been 
proved for WG or equivalently can be easily proved for the VEM (IT8l > on v h ° , we obtain a 
VEM on V® with optimal order of convergence in both H 1 and L 2 norms. 

In other words, we can view a version of VEM as the Schur complement equation of 
WG by eliminating high momentum term and apply the convergence theory developed for 
WG to VEM. 

For elliptic equations with low order terms, e.g., the reaction-diffusion equation, the 
L 2 -inner product of functions in WG is computed as 

{u 0 ,v 0 ) = (n fe (Iit 0 ),n fc (Ir;o)). 

While in VEM, it is computed as 

(nfcM,ngu), 

where I ijj is defined in JT] and can be computed using IljL In WG, however, it is more 
natural to use uo- 

4 . Weak Galerkin Methods and Non-conforming VEM 

In this section we review the original version of weak Galerkin methods 0EI and show 
the equivalence to non-conforming VEM J4). 

Given a polygon mesh 7h and an integer k > 1, we introduce 

w h = {v = {v 0 ,v b },v 0 \K e P k (K),v b \ e G Pfc_i(e) Ve C K,VK c %,}, 

w£ = {vew h ,v b \on = 0}- 

As we mentioned before, now v b is also discontinuous at vertices. The d.o.f. x v can t> e 
dropped and the edge d.o.f. is upgraded to the edge momentum to the order k — 1. The 
interior part vq will be still determined by volume momentum to the order k. We list the 
d.o.f. below 


X = {x 6 ,x 0 } : = {x k e 1 ,Xk}- 
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The weak gradient V w to VPfc is defined as before. The stabilization is changed to 

hft (Ub - QbU 0 ,Vb - QbVo)dK, 

where recall that Qb is the L 2 -projection from L 2 (e) to Pfc_i(e). A more convenient 
formula is 

Xb( u b - Uo\ok) ■ Xb( v b ~ t'olojf)• 

The i 2 -projection Qi, is simply replaced by edge d.o.f. 

The WG finite element is: find Uh G W) j 1 

(19) a WG {u h ,v) = (f,v o) Vv = {v 0 ,v b } G W%. 

where 

« WG K)Vfi) : = V w v h ) + ^ Xb( u b - u 0 \ok) • Xb( v b - v 0 \ dK ). 

KdT h 

The non-conforming VEM space is defined as follows. First the local space defined as 
V k (K) = {v£ H\K) : Vv-n G P fc _i(e),Ve C dK, Av G P fc _ 2 (iT)}. 

To glue them to get a non-conforming space of H 1 , we define 

H 1,nc (Th',k) := (u G H l {K)yi< C Th, f vp ds is continuous Vp G Pk-i(e),Ve G £h}■ 

J e 

The global non-conforming VEM space is 

V h = {v h G H^iTh', k) : v h \ K G V k (K) VK G %}, 

V® = {v h G V h , f vpds = 0 Vp G Pfc-i(e), Ve C dfi}. 

J e 

The d.o.f. of non-conforming VEM space will be the momentum on edges up to order 
k — 1 and the volume momentum up to order k — 2. Again we could enlarge to Vh by 
asking Av G P/. : and including volume momentum to the order k in d.o.f. 

The non-conforming VEM is: find Uh G V)} such that 

(20) «r M K, V h ) = (fh, v h ) Vv h G V; J°, 

where 

ah EM (u,v ) := (Vn^,Vn^)+ 

KeT h 

and for k > 2 

(fh,V h ) = (Qk-2f,Qk-2V h ), 

with a similar treatment for the lowest order case k = 1; see, e.g., cos. 

Just as before, we could recast WG as a VEM in the larger space Vh- The additional 
d.o.f. can be eliminated from the stabilization term. The non-conforming VEM is: find 
Uh G V° such that 

(21) 4 EM K, v h ) = (/, I l k v h ) Vv h G V°, 

where 

ah EM {u,v) := (vn^u,vn^)+ ^ X 6 (m - U k u) ■ x b {v - n fe u). 

KGT h 

The interpolant 11 is defined as before. 


EQUIVALENCE OF WG AND VEM 


11 


Theorem 4.1. The modified non-conforming VEM method (12lb is exactly the original WG 
method (HU. The original non-conforming VEM (120b is the Schur complement ofWG (ED 
by eliminating the high momentum part. 

Advantage of the original formulation of WG and non-conform VEM is the generality in 
aribitary dimension while the conforming version requires the extension of vertices valued 
d.o.f. to polynomials inside faces which now is a lower dimensional polytope. 
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